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Abstract 

The correspondence limit of the atomic elliptic state in three di- 
mensions is discussed in terms of Nelson's stochastic mechanics. In 
previous work we have shown that this approach leads to a limiting 
Nelson diffusion and here we discuss in detail the invariant measure 
for this process and show that it is concentrated on the Kepler ellipse 
in the plane z = 0. We then show that the limiting Nelson diffu- 
sion generator has a spectral gap; thereby proving that in the infinite 
time limit the density for the limiting Nelson diffusion will converge 
to its invariant measure. We also include a summary of the Cheeger 
and Poincare inequalities both of which are used in our proof of the 
existence of the spectral gap. 

1 Introduction 

Nelson's stochastic mechanics [11, 12, 13] encapsulates the dynamical content 
of the Schrodinger equation. It describes the motion of a quantum particle of 
mass m using a diffusion process driven by a Brownian motion with diffusion 
coefficient h/2m. This diffusion process obeys a stochastic generalisation 
of Newton's second law of motion known as the Nelson-Newton law; the 
mass times the stochastic acceleration equals the force. Considering the 
dynamics of such a Nelson diffusion leads to a set of equations which are 
equivalent to the Schrodinger equation. Consequently, this approach allows 



the development of the theory of quantum mechanics without the need to 
entirely abandon the classical world. The fundamental aspects of quantum 
theory such as the uncertainty principle are automatically built into this 
stochastic framework. Moreover, Nelson's theory gives a natural setting in 
which to consider the correspondence limit transition between the quantum 
and classical worlds. 

In [9] we investigated the Nelson diffusion for the atomic elliptic state in 
the Coulomb potential. We derived a limiting Nelson diffusion according to 
Bohr's correspondence principle and showed that in a two dimensional set- 
ting, the trajectories of this limiting diffusion converged to those of a particle 
describing Kepler's ellipse at a rate consistent with Kepler's laws of plane- 
tary motion. This solved a long standing problem in quantum mechanics in 
the two dimensional case. In this paper we will show how these results can 
be extended to the three dimensional case. Our main result is to show that 
there exists a spectral gap for the generator of the limiting Nelson diffusion 
for the atomic elliptic state. It follows that in the infinite time limit this 
process will converge to its invariant measure. Moreover, we will show that 
the invariant measure is, in the limit, concentrated in a neighbourhood of 
the Kepler ellipse in the plane z — 0. 

We begin by recapping the construction of the limiting Nelson diffusion 
from [9]; in Section 2 we give an account of Nelson's stochastic mechanics 
before we introduce the atomic elliptic state in Section 3 and discuss the 
correspondence limit in Section 4. In Section 5 we discuss the invariant 
measure for the limiting Nelson diffusion and show that it is the physically 
correct measure concentrated on the Kepler ellipse in the plane z — 0. In 
Section 6 we discuss Cheeger's inequality and the Poincare inequality which 
we use to prove the existence of the spectral gap in Section 7. 

2 Nelson's stochastic mechanics 

We begin with a brief look at the dynamics of diffusion processes [12] . Con- 
sider a particle of mass m diffusing through M. d according to the Ito stochastic 
differential equation, 

dX(t) = b(X(t),t)dt + edB(t), t > 0, (1) 

where B = (B\, B2, . . . , Bo) is a <i-dimensional Brownian motion with nor- 
malisation E {Bi(t)Bj(s)} = 5ij mm(t, s) for i,j = 1,2, ... ,d. It follows from 
Ito's formula that the density p(x,t) defined for all Borel sets Actf by, 

¥{X(t)eA}= [ p(y,t)dy, 

JyeA 
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satisfies the forward Kolmogorov equation, 

§£ = V- (jVp-bp^=Q*p, P (x,0)= Po (x) (2) 

where Q* is the L 2 adjoint of the generator Q for the diffusion X: 

g= e l A + b . V , (3) 

The paths of X are nowhere differential) le with probability 1, and so to 
discuss the kinematics of diffusions, Nelson introduced the mean forward and 
backward derivatives as conditional expectations, 



D±f(X(t),t) := limE 



f(X(t±h),t±h)-f(X(t),t) 



X(t) 



hlO { ±h 
The mean forward and backward drifts b± are defined by, 

b + (X(t),t) := D + X(t) = b(X(t),t), (4) 

and 

b-(X(t),t) := D-X(t) = b(X(t),t) - e 2 Vlnp(X(t),t). (5) 
The osmotic and current velocities u and v are defined by, 

v:=±(b + + b_). (6) 

The dynamics are introduced through the stochastic acceleration, 

a(X(t),t) = l - (D + D_ + D_D + ) X(t), 

and it follows from Ito's formula that, 

a(X(t),t) = (j£ + (v ■ V)« - (u ■ V)« - j Auj (X(t), t). (7) 

We will now discuss the connection between the dynamics of diffusion 
processes and quantum mechanics. Consider a particle of mass m in the 
quantum state ip = ip(x,t) satisfying the Schrodinger equation, 

ih d $ = -— A^P + Vip, xeR d , t>0. (8) 
at 2m 

We can write this wave function in the form ip = exp(R + iS) where R and S 
are real valued. Nelson's theory states that such a particle can be considered 
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to diffuse through R d according to the Ito stochastic differential equation (1) 
where we choose b = e 2 V(i? + S) and e = ^h/m. We refer to such a process 
as a Nelson diffusion. This can be seen by writing the Schrodinger equation 
(8) in the form, 

l nr^ = -^ i r^ + vm 2 , (9) 

where ip* is the complex conjugate of ip. Equating imaginary parts gives the 
continuity equation for the probability current j, 

|£ + V • j = 0, j := ~ W) = e 2 P VS, (10) 

where p is the particle density p := \ip\ 2 = exp(2i?). The continuity equation 
(10) can be written as, 

^ = V • ^yV exp(2i?) - e 2 exp(2R)V(R + 5) 

which is the forward Kolmogorov equation (2) for the Nelson diffusion with 
drift b = e 2 V(R + S) and density p = exp(2R). Therefore, we see that the 
Nelson theory is consistent with Born's probabilistic interpretation for the 
quantum mechanical particle density p. 

It now follows from equations (4), (5) and (6) that, 

b + = e 2 V{R + S), b_ = e 2 V{S - R), u = e 2 VR, v = e 2 VS. 

Moreover, from equation (7) and the real part of equation (9) we obtain the 
Nelson-Newton law, 

ma(X(t),t) = -VV(X(t),t). (11) 

This is the dynamical content of the Schrodinger equation expressed in terms 
of the Nelson diffusion. 

Since ip satisfies a complex heat equation, it follows that S = —ie 2 In ip 
must satisfy a complex Hamilton Jacobi equation with viscosity ie 2 , 

^! + W + ^ = ^A5. (12) 
dt T 1 m 2 y J 

Hence, defining Z by a complex Hopf-Cole transformation, 

Z = -ie 2 Vln*P = e 2 V{S-iR), (13) 
we obtain the complex Burgers equation, 

f + ,Z.V)Z + ^AZ, (U, 



where the drift of the Nelson diffusion is given by 

b = Re(Z) -Im(Z). 

As is well known, equating real and imaginary parts of the complex Burg- 
ers equation (14) gives the Madelung fluid equation and its continuity equa- 
tion for the current velocity v = e 2 VS, 

dv e 4 /Ap 1/2 \ dp 

- + (v .V)v + VV=-v[^y ^ + V.0») = 0. (15) 

When if; is a stationary state we can write, 

ip{x,t) = iPe{&) exp(—iEt/h), t > 0, 

and then i/je(x) = exp(R + iS) which satisfies the stationary Schrodinger 
equation, 

h 2 

-—AtJj e + V4> e = E^ e . (16) 

Then Z satisfies, 

1 , V it 2 

_ z2+ ___ V .Z = £ . (17) 

Moreover, for a stationary state, p = exp(2R) is the invariant density of the 
quantum system, and so is also an invariant measure for the Nelson diffusion. 



3 The atomic elliptic state 

We begin with the atomic circular state |Circ) n for a particle of unit mass 
given by, 

(cc|Circ) n = *„,„-i,„-i(aj), 

where x = (x, y, z) and ^ n ,i,m (with I — 0, 1, . . . , n— 1, \m\ < I and m, n G N) 
is the nodal Schrodinger wave function for the Hamiltonian of the hydrogen 
atom, 

with [qk,Pi] = ihSki for k,l = 1,2,3. We choose suitable units so that h = 1 
and /i = 1. Then, for the orbital angular momentum L, where 

L = (Li, L 2 , L 3 ) = q X p, L 2 = (L\ + L\ + Lf), 
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we have the eigenvalue relations, 

The atomic elliptic state is then given by, 

|EUiptic(0)) n = exp(-iM 2 )|Circ) n , 

where, 

A = (A U A 2 ,A 3) - 1 ^Xi-rxp) , 



v /= 2E V 2 l«l 

is the Hamilton-Lenz-Runge vector. 

Following [10], to within a multiplicative constant (after reinstating h and 
the force constant //), 

i/> e , n {x) := (x | Elliptic (0)) n = exp — |a:|J £ n _i(ra/), - (18) 

where 

with as = (x,y,z), E n = — /i 2 /(2A 2 ), e 2 = h, A = ne 2 and £ n _i a Laguerre 
polynomial. This wave function is localised on the ellipse with eccentricity 
e = | sin 6* | given by, 

where a = A 2 //z. We call £ e the Kepler ellipse. For integer n this wave 
function satisfies the stationary Schrodinger equation, 

-^e 4 A^ n (x) - ^e,n(x) = E n ^ n {x). (20) 

Remark 1 . Throughout this paper e will refer to the eccentricity of the Kepler 
ellipse. When we use the exponential function we shall always write exp(x) 
rather than e x to avoid confusion. 
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4 The correspondence limit 

4.1 The limiting state and limiting Nelson diffusion 

Let the atomic elliptic state from equation (18) be written in the form, 

ipe : n(x) = exp(i2 ej „ + iS £:n ). 

This stationary state has a Nelson diffusion X e , n (t) given by, 

dX etn (t) = b e>n (X e>n (t)) dt + e dB(t), 

where b e ^ n = e 2 V(R^ n + S €jn ). For integer n, the diffusion X e ^ n satisfies the 
Nelson-Newton law (11). 

We now take the correspondence limit of X e;n by considering the drift 
b etn as n — > oo while e — > with A = ne 2 fixed. This gives a new diffusion 
X e such that, 

dX e (*) = &(X e (i))df + edB(i), 



where 



b= }im e 2 V(R £ , n + S e , n ). 

n|oo,e4.0 
A=ne 2 



We will refer to X e as the limiting Nelson diffusion. 

For this we use the vector Z e>n = — ie 2 V In ^ £jTl introduced in equation 
(13). In Cartesians, 

A \ L n -\{nv) ) \x\ Xe C n -i{nu) \ J 
where v is as defined in equation (19). 
Lemma 1. Define, 



Zq,oo '■— hm Z en . 

nT°o,eJ.O ' 



Then, 



where x = (x,y, z) and v is as defined in equation (19). 

Proof. See [9] □ 
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Lemma 2. Let, 
Then, 

VR e = 



Z Qt00 {x) = e 2 V{S t -iR e ). 



^{(l + a)^ + ((l-a) >/3 ^0)} 

vs. = -^(^(-Ad-^o)}, 



where 



a 



2 





X 




X 


- 4A 2 e//x) 2 + (1 - e 2 )y 2 


(e| 


X 


- x) 2 + (1 - e 2 )y 2 



+ 2 







- x - 2A 2 e/^) 2 + (1 - e 2 )y 2 - 4A 4 e 2 //i 2 


(e| 




- x) 2 + (1 - e 2 )y 2 



p 



-2\ 2 eVT^y 



\x ((e|jc| — re) 2 + (1 — e 2 )y 2 ) a' 
Proof. A simple calculation from Lemma 1. 
By setting, 

b = lim e 2 V(# £ , n + S e , n ) = Re(Zo i00 ) - Im(Z 0iOO ) = e 2 V(# e + 5 e ), 

n|oo,ej.0 

we find the following limiting Nelson diffusion: 

Theorem 1. The limiting Nelson diffusion X t satisfies the ltd equation, 

dX e (s) = b(X e (s))ds + edB(s), 
where b(x) = (b x ,b y ,b z ) in cartesian coordinates with, 

K = A{(„ + /> -i)i-(„ + /> + i)* 



6 » = 2x{ (a " /3 " 1) ^7^"" (a + /J+1) jlj 

with a and (3 as defined in Lemma 2. 



Proof. See [9]. □ 

As a result of taking this limit, X e does not satisfy the Nelson-Newton 
law and so is not a true Nelson diffusion. This is because it is not associated 
with an exact solution to the stationary Schrodinger equation. However, if 
we let ip t = exp(_R e + iS € ), we find an approximate wave function which is 
the correspondence limit of the atomic elliptic state. 

Theorem 2. The correspondence limit of the wave function ip^ n is the lim- 
iting wave function ip e , where 



v 




and v is as defined in equation (19). 

Proof. See [9]. □ 

Simulations of the limiting Nelson diffusion are shown in Figure 1. The 
paths all rapidly converge to a neighbourhood of the Kepler ellipse. 

Remark 2. As expected from equation (17), the vector Z 0>oo satisfies, 

2 Z °'°° 1*1 " 2A 2 ' 1 ] 

Consequently, by taking imaginary parts of equation (22) we find, 

VR e ■ VS e = 0, (23) 
which will be key in what follows. 




Figure 1: Simulations of the limiting Nelson diffusion X e for e = 0.5, e = 0.1. 
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4.2 The drift singularity 



The limiting Nelson diffusion X e has a singularity in the drift field b which 
arises from the existence of nodes in the atomic elliptic state wave function 
ip e>n . The nodes of ip e>n are given by the point |ac| = and also by the right 
branch of a series of n — 1 hyperbolas in the plane y = 0. The drift b e n of the 
Nelson diffusion X t n has infinite blowups on these curves. This behaviour 
makes the nodes inaccessible to the Nelson diffusion [2]. The drift b for the 
limiting Nelson diffusion X 6 also has an infinite blowup at the point |a?| = 0. 
However, the correspondence limit of the hyperbolas forms a finite jump 
discontinuity in b across the surface a = 0. This surface of discontinuity is, 



(x,0,z) 



1 - e 2 



< x < 



1 -e 2 



This limiting singularity and the nodal curves for the exact case are shown 
in Figure 2. They are both ringed by the Kepler ellipse. 





Figure 2: The nodal surfaces for ip etJl for e = 0.5 and n — 10 together with 
their correspondence limit £ shown with the Kepler ellipse. 



4.3 Unifying the two and three dimensional analyses 

We will now prove that our results in [9] for the two dimensional case are 
relevant for the full three dimensional problem. Recall that for integer n 
the three dimensional atomic elliptic state ip e>n (x,t) satisfies the stationary 
Schrodinger equation, 

-^e 4 A^ in (aj) - T^Tip ejn (x) = E n if) ei n(x), 
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where e 2 = H, E n = —fj 2 /(2X 2 ), X = ne 2 and x = (x,y,z). 
It can be shown that the two dimensional function, 

«<«> = ex P (-^^h) k* (v^nyf ) , 

where x = (x,y), v is as in equation (19) and Ti2n is a Hermite polynomial, 
satisfies the stationary Schrodinger equation, 

where again e 2 = H, E n = — /i 2 /(2A 2 ), but now A = ne 2 + e 2 /2. 

We will investigate the connection between the correspondence limits of 
these two and three dimensional wave functions. For this we need the fol- 
lowing result on Hermite polynomials. We follow the conventions in [3] . 

Lemma 3. Let, 

i n' m (v^+iv^) 

where Ti m denotes the m th Hermite polynomial. Then, 



ton «2 3= = l- A /l_i. 



A = ne 2 +e 2 /2 

Proof. Clearly, 

Qm / = X ^ ■ \ • 

Then using the standard recurrence relation, 

H m +i{v) = 2vH m {v) - 2mH m -i{v), 
with v = \/m + \\Jvj2 we can conclude, 

Vm + 2 y m+ i 
Assuming that Qm — ► 9 as m — > oo then, 

^ = 2v^72-g, 

giving the desired result. □ 
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Theorem 3. Let ipf d (x,y, z) denote the correspondence limit (in the sense 
of Theorem 2) of the three dimensional atomic elliptic state i/j €tn (x, y, z) and 
let ip 2d (x,y) denote the correspondence limit of the two dimensional wave 
function ijj 2d n (x, y) . Then, 



tf d (x,y)=^ d (x,y,z)\ z=0 . 



Proof. If we define, 
then in Cartesians, 

z 2d Jx) = if I i 



Z 2 t d n {x) = -te 2 Vln^ 2d 




where Q m is as defined in Lemma 3. Taking the correspondence limit gives 



njoo,e|0 
A=ne 2 + e 2 /2 




where x = (x,y). Therefore, from Lemma 1 we can conclude, 

(Zo%(x,y),0) = Z ^(x,y,z)\ z=0 , 
which gives the result. □ 

In [9] we worked with the limiting Nelson diffusion restricted to the plane 
z = which we can now view as being derived from the correspondence limit 
of a two dimensional wave function. We showed that as e — > this two di- 
mensional limiting Nelson diffusion converged to the deterministic dynamical 
system X , where 

X = b(X ), 

provided the path avoided the drift singularity S. Moreover, we showed that 
the Kepler ellipse was an asymptotically stable periodic orbit for this dynam- 
ical system whose domain of attraction contained all of IR 2 outside the Kepler 
ellipse. We also showed that a particle on the Kepler ellipse would move ac- 
cording to Kepler's laws of planetary motion. Thus, we could conclude that 
in the infinite time limit, any trajectory beginning outside the Kepler ellipse 
would converge to Keplerian motion on the Kepler ellipse. These results were 
proved using the following non-orthogonal coordinate system which we will 
again use: 
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Definition 1. The Keplerian elliptic coordinates (u,v) in the plane (x,y) 
are defined by, 



2ae(cost> — u) 2ae\/l — u 2 sin v 

x= , y= ■ , 

e + u e + u 

where — e < u < 1 and < v < 2n with a = \ 2 /fi. 

The u coordinate curves form a family of non-confocal ellipses with eccen- 
tricity foci at (0, 0) and (— Aaeu/(l+e), 0), and semimajor axis 2ae/ (e+u) 
parallel to the x axis. The Kepler ellipse is given by the coordinate curve 
u — e, the singularity S by the curve u — 1 and the ellipse at infinity by 
u = —e. On the Kepler ellipse v is the eccentric angle. 

Unfortunately we have been unable to find an extension for this system 
to simplify the three dimensional limiting Nelson diffusion. However, the 
following naive extension will be of use in studying the invariant measure. 

Definition 2. The cylindrical Keplerian elliptic coordinates (u,v, z) in three 
space are defined by, 

2ae(cosf — u) 2ae\ / l — v? sin v 



where — e < w < 1, < u < 27r, z eR with a = A 2 //i. 

It follows from Theorem 3 that if we could prove that the motion was 
confined to a neighbourhood of the plane z = then our two dimensional 
results would be relevant. This provides the motivation for studying the 
invariant measure for the full three dimensional problem. 



5 The invariant measure for X e 

We again recall that the atomic elliptic state ip^ n = exp(R en + iS e , n ) satisfies 
the stationary Schrodinger equation (20) provided n is an integer. There- 
fore, assuming n is an integer, the Nelson diffusion X e>n satisfies the Nelson- 
Newton law and has invariant density p e ,ri = exp(2i? ein ). Thus the forward 
Kolmogorov equation (2) gives, 

= AS e , n + 2VR e , n ■ VS e , n . 

Moreover, from equation (18) it is possible to write, 
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where R e>n and S ttH only depend on e through A = ne 2 . Thus, 

= e 2 AS e , n + 2VR e , n ■ VS e , n . (24) 

This equation is equivalent to p e „ being the density for the invariant measure 
for X e>n . 

The limiting Nelson diffusion X € does not satisfy the Nelson-Newton 
law and the limiting state ip e = exp(i? e + iS e ) is only an approximate wave 
function. Therefore we cannot automatically conclude that exp(2i? e ) is the 
invariant density for X t . However, we can write R t = R/e 2 and S € = S/e 2 
and then from equation (24) conclude that exp(2i? e ) would be the invariant 
measure if and only if, 

= e 2 AS + 2VR ■ VS. (25) 

Thus, from equation (23) with the observation AS ^ we find that exp(2i? e ) 
is the invariant measure for the limiting Nelson diffusion only when e — > 0. 
We now show that this is the physically correct invariant measure. 

Lemma 4. Let the three dimensional limiting wave function be given by 
ip e = exp (R e + iS e ) where R e and S e are real. Then, 

OR 

-^(x,y,z)=0 & z = 0. 

Proof. Follows from Lemma 2. □ 

Theorem 4. R e (x,y,z) has a manifold of constant maxima on the Kepler 
ellipse in the plane z — 0. 

Proof. From the previous lemma we know that z = is a necessary condition 
for R £ to have a critical point and so we can write a and (3 from Lemma 2 
in Keplerian elliptic coordinates as, 

a/(1 — u 2 )(l — e 2 ) (e + u) sin v 

OL = ~ } \ , P ~ 



1 + eu — (e + u) cos v J 1 + eu — (e + u) cos v 

The Kepler ellipse is the curve u = e and so again from Lemma 2 we find 
that V-R e | 2=0u=e = 0. Moreover, it can easily be shown that on the Kepler 
ellipse, 

so that R e is constant around the ellipse. Clearly, any point on the Kepler 
ellipse is therefore a degenerate critical point of R e forcing the Hessian of R e 
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to be singular. However, in cylindrical Keplerian elliptic coordinates (it, v, z) 
the Hessian on the Kepler ellipse is, 



K(%v,z)\ z=0iV=e 



/ A(l+e 2 +2ecos?j) 
4e 2 e 2 (l-e 2 ) 





V 


















e 2 A 3 (l- 



-2ecosw) 



/ 



confining the degeneracy to the v coordinate as expected. It is now clear 
that this must be a manifold of maxima. □ 

This maximum forms a volcano shape in the plane z = for exp(2i? e ) 
as is shown in Figure 3. By considering the asymptotic behaviour of this 
function as e ~ we can find the thickness of this volcano in the z = plane 
and also the thickness in the z direction: 




Figure 3: exp(2i? e ) for z = 0.2, z = 0.1 and z = with e = 0.5 and e = 0.1. 



Theorem 5. For an arbitrary test function G G ^(1R 3 ), 
as e ~ 0, where 

g(v) = (1 — e cos v ) y/l + e 4 — 2e 2 cos 2v, 

G(v) = G(x(v)), x(v) = (a cos v — ae, ayi — e 2 sint>, 0), 
and v is the eccentric angle of x(v) G S e . 
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Proof. We can write the integral using cylindrical Keplerian elliptic coordi- 
nates and then apply Laplace's method twice, 



G(x) exp 



JO J-e J - 



2R(x) 



dx 



G(u, v, z) exp 



2R(u, v, z) 



d{x,y,z) 



d(u, v, z) 



dz du dv 



G(e, v, 0) exp 1 



d(x,y,z) 



d(u, v, z) 

The result then follows from a simple calculation. 
Corollary 1. As e ~ ; 



Tie 2 


/ Tie 2 


d 2 R\ \ 


/ \d 2 R 


dz 2 \ \ 


1 \ du 2 



dv. 



z=0,u=e 



□ 




exp | ] ,\ x 



16 \^ 27ra 3 e 2 



J \e* J A 

where E(£) denotes the complete elliptic integral of the second kind and 

Proof. Follows directly from Theorem 5. □ 
Corollary 2. The effective standard deviation of the Gauss measure, 

exp 1 2R ^> 



measured normal to the Kepler ellipse in the plane z = is, to leading order 
in e, 

e\z I (1 — e cos v) (1 + e 2 + 2e cos t> ) 



and m t/ie z direction is, 



1 + e cos v 



e\2 



2e cosf . 



Proof. Combining Corollary 1 with Taylor's formula gives, 



d 2 R 



{x e )hihj + 0(e) 



(26) 
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for some constant k where x e is a point on the Kepler ellipse and h = 
(hi, hi-, h?) is a vector in the plane through x e orthogonal to the Kepler 
ellipse. Since R has a maximum on the Kepler ellipse, its Hessian is negative 
definite and so the leading term in equation (26) can be written as, 



cxp 



\h\ 



E 



d 2 R 

dxidxj 



(x e )hih 3 



where h is the unit vector parallel to h. Therefore, we can conclude that the 
width of the volcano in the direction h can be measured in units, 



E 



ho 



d 2 R 

dxidxj 



(x e )hih 3 



giving the result. 



□ 



We now look to find an invariant density for X t when e 7^ 0. We follow 
the methods of Truman and Zhao [15] (for zero noise) and make the ansatz 
that the invariant density can be written in the form, 



p e (x) := T(x) exp 



2R(x) 



cxp 



2R(x) 



+ \nT(x] 



(27) 



where T is 0(1) in e. As e ~ 0, the exponential term will vary rapidly in 
space variables, with a pronounced peak of constant maxima on the Kepler 
ellipse in the plane z — 0, while T will only vary slowly in space variables, this 
variation being in the direction tangential to the Kepler ellipse. Let n t denote 
the stationary probability measure with density (before normalisation) given 
by p e . We can show that this invariant measure is concentrated on the Kepler 
ellipse for small e: 

Corollary 3. For an arbitrary test function f G y(M. 3 ), i/E denotes expec- 
tation with respect to ir e , 



Cf(v)T(v)g(v)dv 



fo* T ( v )9(v)dv 



as e ~ 0, where f(v) = f(x(v)) with x(v) G £ e and g(v) as in Theorem 5. 
Proof. Follows from Theorem 5. □ 
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We can also find the function T up to leading order in e 2 . Assuming p e 
is the invariant density, the forward Kolmogorov equation (2) yields, 



6 2 



-VS+-VlnT)p e ) =0, 



giving, 



e 2 

-VlnT- (VS- Vi?-e 2 VlnT) - AS + -A InT = 0, 

which from equations (5) and (6) can be written, 

VlnT b_ +V v = 0. 

If we take only the leading term in e 2 then, 

{VlnT)-V(S -R)+AS = 0. 

Assuming the variation in T is parallel to the Kepler ellipse we can write 
V In T = 7b as the drift b is tangential to the Kepler ellipse when the particle 
is on the ellipse. Moreover, on the Kepler ellipse V-R = and so, 

VS 

VlnT = -AS-, 



|VS| 2 ' 

Using the Keplerian elliptic coordinates on the Kepler ellipse so that v be- 
comes the eccentric angle, we see that the variation in T as the particle 
precesses around the Kepler ellipse is, 

d -lnT = VlnT dxM ^ ^ 



dv ' dv 1 + e 4 — 2e 2 cos 2v ' 

Hence, 

. ^ ( C 2e 2 sin2£ , \ (1 - e 2 ) 

T(v) = exp - / — — — dv * - 



o l + e 4 -2e 2 cos2{} J Vl + e 4 - 2e 2 cos2t;' 

This gives the slowly varying multiplicative behaviour of the invariant density 
p e which is shown in Figure 4. This behaves exactly as one would expect from 
the results in [10]. 

Corollary 4. For an arbitrary test function f G y(M. 3 ), i/E denotes expec- 
tation with respect to n e , 

E{/} ~2^/ f(v)(l- ecosv) dv, 
as e ~ 0, where f{y) = f(x(v)) and x(v) G £ e . 
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Proof. Follows from Corollary 3 with, 

T(v)g(v) = (1 — e 2 )(l — ecosf). 



□ 




Figure 4: The function T(v ) with the invariant density p e . 



6 Some preliminary results 

Before we discuss the spectral gap for the generator of the limiting Nelson 
diffusion process, we will first introduce some of the tools for our proof; the 
co-area formula, Cheeger's inequality and the Poincare inequality. The first 
two of these are standard results in Riemannian geometry which we will 
recast in Euclidean space. We let M C M 3 be an open Borel set and let \i be 
a probability measure on the Borel sets of M with a smooth density p with 
respect to 3-dimensional Lebesgue measure. Let po be the 2-dimensional 
measure with density p with respect to 2-dimensional Riemannian volume 
on surfaces. 

Lemma 5 (An elementary co-area formula). Let f : M — > R be a smooth 
measurable function. Then , if m := sup xgA/ f 2 (x), 




Proof. Consider the level surfaces of the function f 2 (see Figure 5). Let 
n = V//|V/| denote the unit vector orthogonal to the surface f 2 = s. Then 
for small o > 0, 

f 2 (x + an) = s + ds, 
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Figure 5: The elementary co-area formula 
and so to first order, 

V(/ 2 ) • an = ds. 

Evidently d/i — a dfig giving, 

dfi d ds = 2dfi d a\f\\Vf\ =2d/i|/||V/|, 

and the result follows. □ 

The following inequality due to Cheeger [6] originates in the study of 
lower bounds for the eigenvalues of the Laplacian on a compact Riemannian 
manifold with Dirichlet boundary conditions. Our sketch proof is based on 
the presentations of Buser and Chavel [4, 5]. 

Theorem 6 (Cheeger's inequality). Let, 

k = mf 

A pL(A) 

where A C M varies over all open sets with smooth boundary satisfying 
A fl dM = 0, and let f : M —> M. be a smooth measurable function with a 
suitable boundary condition e.g f\dM = 0. Then, 

(MM>V 

(//'d/x) -T- 

Proof. Let m := max^gM f 2 {x), then by the co-area formula (Lemma 5), 

pm pm 

Vf\dfi= fi 9 (f 2 = s)ds>h t,(f 2 >s)ds, 

M JO JO 
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where 

h : = mf — — 

0<s<m fi[f 2 > S) 

Integration by parts gives, 



and, 
Thus, 



M/ 2 > 8 )ds=-J (/i(/ 2 > 8)) ds, 



f s (M/ 2 > s)) = -Uf = s). 



M 



\Vf 2 \dfi>h I pi(f 2 >s)ds = hj spL 9 (f 2 = s)ds 

= h I fd^ d . 



The Cauchy-Schwartz inequality gives, 



j fd^ (| iwiM/i) > (| /iv/id^y, 



and so, 



(/ | V/| 2 d/x) (J/lV/ld/i) 2 = (JlV/^ld/x) 2 h* 



Finally, it is clear that h > h giving the result. □ 

Our proof also uses the celebrated Poincare inequality and its perturba- 
tion theory. In particular we use the following classical Poincare inequality 
for the Lebesgue measure of open balls in R 3 . 

Theorem 7. Let dx denote integration with respect to the 3 -dimensional 
Lebesgue measure. For all f G C 1 ( 



fdx<% [ iv/fda+^/y fdx)\ 

B r ^ JB r \ B r\ \JB r J 

where \B r \ denotes the Lebesgue measure of the ball B r . 

Proof. A proof of the Poincare inequality on Euclidean balls can be found in 
[14]. The constant 4/7T 2 used in the inequality is derived from Theorem 1.3 
of [16]. □ 

We are now ready to prove the existence of the spectral gap for the 
generator of the limiting Nelson diffusion. 
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7 The spectral gap 



Consider any Nelson diffusion process X with drift b. The generator of this 
process is, 

g= € -A + b-V. 

We can rewrite this operator using a simple similarity transformation. 
Lemma 6. For the generator Q of any Nelson diffusion process, 

G = -^exp(-i?- S)Hexp(R + S), 

where H is the formal Hamiltonian, 

H=l(-e i A + e 2 V-b + b 2 ). 

Moreover for ip = exp(i? — S), 

Htp = 0. 

Proof. Simply follows from repeated application of the identity, 
A(fg) = fAg + gAf + 2Vf-Vg, 

together with the condition b = e 2 V(i? + S). □ 

Unfortunately, we have not been able to prove for the generator Q of the 
limiting Nelson diffusion X e defined by, 

g = ^A + b-v, 

where b is as in Theorem 1, that H has a spectral gap as is explained in 
[9]. Therefore, we have to resort to different techniques to understand the 
spectrum of this generator. Our aim here is to establish a suitable Poincare 
inequality which will give us the existence of the spectral gap. This inequality 
is derived from separate inequalities which are established both inside and 
outside some open ball following a decomposition argument as in [17]. A 
detailed discussion of the use of Poincare inequalities with diffusions can be 
found in [18]. 

For any r > 0, define the open ball B r = {x E 1R 3 : \x\ < r} and 
its complement B c r = {x e IR 3 : \x\ > r}. Let Jzf (-Ev) denote the set of 
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all Lipschitz continuous functions with support in B^. We assume that the 
invariant density for the limiting diffusion X e can be written as, 

p e = T{x)exp{2R e {x)) 

and that there exists a constant r > such that, 

|VlnT(aj)| < C (28) 

for all x with \x\ > r where e G (0, 1) is sufficiently small so that, 

0<C<£. 

We also assume that In T(x) is bounded for all x e M 3 . If we define V e = 
2R e + \nT then we have dir e = exp(V e ) dx. 
Define a new operator, 

0„:=yA + (t*.V), 

where u = (e 2 /2)Vlnp e is the osmotic velocity in the stationary state corre- 
sponding to the invariant measure with density p e . 

Lemma 7. For any Borel set A C M 3 , 

- [ g u \x\d7r e = ^- [ n- \/\x\d(7r t ) d , 

J A 1 JdA 

where n is the inward pointing unit normal to the surface dA. 
Proof. Since, 

p e A|aj| — \x\Ap e = V • (p e V|cc| — |cc|Vp e ) , 

and, 

V • (\x\p e u) = \x\V ■ (p e u) + p e u ■ V\x\, 

it follows that, 

fe 2 \ e 2 

peQ u \x\ = \x\ ( — Ap e - V • (p e u) J + —V • (p e V|x|). 

Moreover, from equation (15) it follows that = V ■ (p £ v), where v is the 
current velocity in the stationary state. Therefore, 

e 2 . _ , .A e 2 _ , _, „ e 2 



PeQ u \x\ = \x\ (^-A Pe - V • ( Pe b)j + -V • (p e V\x\) = -V • (p e V\x\), 
and the result follows by Green's theorem. □ 
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Lemma 8. There exists a constant r\ > such that for all x with \x\ > T\, 

Q u \x\ < — — < o, 

where C — (/i — e 2 AC)/(e 2 A) and C is as defined in equation (28). 
Proof. Clearly, 

e 2 

g u \x\ = — (2 + 2VR e -x + VlnT-a;). 
2 \x\ 



Moreover, from Lemma 2, 



and so, 



lim a — 1, lim (3 = 0, 

|cc|^oo |cc|^oo 



lim ^(i + V i2 e . aJ ) = -f. 

|a;|^oo \X\ A 



Therefore, there exists r\ > r such that for all x with \x\ > n, 



i^-r (1 + Vi? e • x) < 
\x\ 



2A' 

which with the assumption in equation (28) gives the result. □ 
Lemma 9. For every g e ^(B^), 



Jg 2 d7r e <i- 2 J \Vg\ 2 dn e 



Proof. Combining Lemmas 7 and 8, 

€^C f €^ f 6^ 

V^( A ) ^ - / d ^ = T / n • V l a; l d (^)^ ^ ^(^)a(^), (29) 

where n is the inward unit normal vector field of dA. It follows from equation 
(29) that for all A C B° rv 

W#>C>0. (30) 

Therefore, using Cheeger's inequality (Theorem 6) on M = B^ with a Dirich- 
let boundary condition on dM = {x : \x\ — r±} with equation (30) gives for 
sufficiently smooth g : B c ri — > K. such that p(rr) = for |x| = r l5 

(/^ 2 d7r e ) " 4 " 4 ' 
giving the result for all <? G Jz^-EvJ. □ 
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Lemma 10. There exists a constant r 2 independent of e with r 2 > r\ such 
that, 



Proof. It follows from equation (30) that, 



- Tr ky(n e (B r ))- ^ m >C, 

and so integrating, 

7r £ (B c r ) < tt^J exp (-C(r - n)) < exp (-C(r - n)) , 

gives the result. □ 
Theorem 8. For any f G Cg(R d ), 

where r% = r 2 + 8/C. 

Proof. For some constant 5 > 0, define, 

{0, |x| < r 2 , 

/(a-), r 2 < |cc| <r 2 + <5, 
/(«), r 2 + <5<|a3|, 

where r 2 is as in Lemma 10. 

It follows that as r 2 > ri then g £ Sf (B^) and so from Lemma 9, 

L f dne = L 92 d7Te - \ 92d ^-h\ 1 v ^' 2 dne - 

Moreover, from the definition of g, 

0<|Vs| 2 <2|V/| 2 + |/ 2 , 

and so, 

/ /^,^/|V/|^, + ^//^. 

r-2+o 

Finally choosing S — 8/C gives the result. □ 
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We now prove a similar result on the interior of the ball using the Poincare 
inequality on Euclidean balls. 

Theorem 9. For all f e C\R d ) with J f dn e = 0, 

J b fd^^exp^ + l)) J |V/| 2 d7r e + | J f dn e . 

Proof. Recall V e = 2R e + InT so by Lemma 2, there exists a constant k > 
such that for any r > 0, 

\V e (x)\<^(\x\ + l), x E B r , ee(0,l]. (31) 



Clearly 



and since dn e = exp(K) dec, it follows from equation (31) and Theorem 7 
that, 



< ^T ex P(^ + 1 )) / |V/| a ds 







exp 


V2e 2 


4r 2 






exp ^ 










7T 2 


exp f 



Therefore, 



It follows from J /d7r e = that, 



2 



(32) 



2 / \ 2 



/d7T £ J = ^ /d7T e 

and so by the Cauchy-Schwartz inequality and Lemma 10, 



2 

p2 . 



/d7T e < / d7T e / /M7T e < - / / d7T e . 



S, 3 / V^B= 3 / V^S? 3 
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Therefore, by (32) we obtain, 

J b / 2 d7T e <^|exp(^(r3 + l)) J |V/| 2 d7T e + | y J 2 d7T e , 

which gives the result. □ 

We can now combine these results to give a single Poincare inequality. 

Corollary 5. For any f G Co(lR d ) such that J f dn e = there exists a 
constant k, such that, 



/ 2 cbT e <exp (k/e 2 ) J \Vf\ 2 dn e 



Proof. Follows directly from Theorem 8 and Theorem 9. □ 

All that remains now is to show that this Poincare inequality leads to the 
spectral gap. 

Theorem 10. There exists a constant 7 = 7(A, fi,e,C) > such that the 
generator of the limiting diffusion process X e , 

£:=^A + 6-V, 

has a spectral gap where, 

gap(£) > exp (-7/e 2 ) , 

for any e G (0, 1]. 

Proof. Let Ai = inf(er(— Q) \{0}) where cr(-) denotes the spectrum of a linear 
operator. Then there exists a sequence of functions {f n }n>i such that, 

//nd7T £ = 0, Jf 2 dlT e = l, 

and, 

lim [(\lfn + Gfn) 2 d7T e = 0. 

By the Cauchy-Schwartz inequality, 

J (Ai/n + Gfn) 2 dn e = Q f n dn^j (J (Ax/ n + Qf n f d7r ( 



> 



J{\ifl + fnSfn) dn^j >0. 
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Therefore, 

Ai = - lim / fnGfnd-Ke. 

Now, for / G V(Q), 



e 2 



fGfdn e = -J f^-Af + b-Vfjp e dx 

= ~\j (^Af-eW + b-V/ 2 )^* 
= ~\J PeQfdx + ^ J \Vf\ 2 dn e 
= -\f f 2 G* Pe dx + ^ J \Vf\ 2 dn e . 
Therefore, as p t is the density for the invariant measure, 

-Jf^fdn^^J IW| 2 dvr e . 

It follows that, 



e 2 



lim — / |V/ n | dn t = Ai. 

n— >oo Z J 

However, from Corollary 5, 

1 = lim / / 2 dvr e < exp(A;/e 2 ) lim / |V/ n | 2 d7r e = ^exp(A;/e 2 ). 
Thus Ai > exp(— 7/e 2 ) for some constant 7 > and all e G (0, 1]. □ 

8 Conclusion 

We have shown that in the three dimensional case for small e > 0, the limiting 
Nelson diffusion derived from the correspondence limit of the atomic elliptic 
state has an invariant measure which is concentrated in a neighbourhood of 
the Kepler ellipse in the plane z — 0. We have also shown that the density of 
the limiting diffusion will converge in the infinite time limit to this invariant 
measure irrespective of the eccentricity of the Kepler ellipse. We can therefore 
conclude that in the infinite time limit the motion of a diffusing particle will 
be confined to a neighbourhood of the z = plane. Hence the results of [9] 
become highly relevant in establishing that this motion becomes Keplerian 
motion on the ellipse in the infinite time limit. We hope to investigate the 
ramifications of this result for planetesimal diffusions [1, 7, 8] in a future 
publication. 
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